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Abstract 

I have investigated what coordinate systems, other than the Robertson- 
Walker (RW) system, there might be to describe a homogeneous and 
isotropic universe that are physical close to the origin, i.e., covariant 
transforms from RW that there have a Minkowski metric. I show that 
all possible such systems must have the same variable photon veloc- 
ity close to their origin, which is consistent with the RW system also 
having the same variable photon velocity. Implausible though it may 
seem, if the criteria for physicality arc valid, it seems inescapable that 
this must represent the physical photon velocity c(t). The assumption 
of homogeneity makes it applicable to physical processes throughout 
the universe. I find that the analytic expression for it is c(t) = ah, 
the cosmic scale factor times the Hubble ratio in appropriate units. 
The existence of variability derives from the physicality criteria being 
applied to the transformed coordinates expanded from the origin out 
to the lowest order of the RW radial coordinate, and not to gravita- 
tional energy density, which determines only its magnitude. I show 
that there are covariant transforms which are physical over a large 
fraction of the universe's space-time. The replacement of the RW time 
by a generalized time, for which the differential is c(t)dt, keeps in- 
variant the Lorentz transform, the Minkowski metric and the Einstein 
field equation. These can then be used to preserve and extend all our 
physical laws. Other physical "constants" that depend on the photon 
velocity may also change. 
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1 Introduction 



The cosmological assumption that the universe is homogeneous and isotropic 
in space does not automatically lead to a constant photon velocity. In fact, 
the Robertson- Walker (RW) space-time line element [1] [2], which can be de- 
rived from the assumption of this symmetry, allows for a variable photon 
velocity c(t) which depends on the cosmic time t. Although there is strong 
historic precedent for assuming a constant photon velocity, it should be no 
more a surprise to us that it is variable than that the universe size is variable 
in time, a great surprise to Einstein when he discovered his equations led to 
it, which he tried to fix with his famous cosmological constant. In this paper 
I propose to show that the same assumptions lead to a variable physical pho- 
ton velocity. I do this by showing that every possible spherically symmetric 
physical coordinate system requires the same variable photon velocity c(t) 
close to the origin. By physical I mean that it is obtained by a covariant 
transform from the RW coordinates and that its coordinates expanded from 
the origin to the first order of RW distance have a generalized Minkowski 
metric, i.e., a Minkowski metric in which c(t)dt replaces codt. If this is a 
valid description of physicality, one would have to conclude that this must 
describe the physical photon velocity. Homogeneity requires each coordinate 
system to be the same no matter where its origin is located in the universe. 
This means that the same variable photon velocity applies to all physical 
phenomena throughout the universe. 

A variable photon velocity has far reaching implications. A large num- 
ber of other physical "constants" depend on it, so many of them presumably 
also change with cosmic time. In addition, a variable c(t) no longer keeps 
the Lorentz transform invariant, so that the equivalency principle must be 
modified. Following the suggestion of Magueijo [3], we introduce a general- 
ized time i for which di = c(t)dt. With this substituted for cot, the Lorentz 
transform, the Minkowski metric, and the Einstein field equations all re- 
main the same for a variable c(t). We will label them "generalized" when 
using t. Of course, this generalized Lorentz (L ) transform and generalized 
Minkowski (M ) metric will reduce to the normal transform and metric for 
events whose time separation is small compared to cosmic time in order that 
all the laws of atomic, nuclear, and particle physics remain valid. Because of 
this, we use this M metric close to the origin as the physicality requirement 
for the transformed metric. 

We should point out that this physicality requirement excludes the proper 
coordinate system. Since the time of the proper coordinates is identical to 
the time of the RW coordinates, there is not the time dilation required by 
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the L transform even though the points on the two systems move at an in- 
creasing velocity with respect to each other as the distance from the origin 
increases. We will require physical transforms near the origin to have at 
least a first order time dilation caused by the L transform. 

We proceed by enumerating in Section 2 the four assumptions needed 
to derive c(t) and three more assumptions to adapt this c(t) to our physi- 
cal laws. We then show in Section 3 that the assumption of isotropy and 
homogeneity leads to RW coordinates which can be generalized to include 
c = c(t). In Section 4 we introduce a general metric for isotropic coordi- 
nates. These are made a covariant transform by finding the 4-velocity of a 
point on the frame of the transformed coordinate system in terms of both 
the RW coordinates and the transformed coordinates themselves. In Sect 5 
we justify using the M metric as the criterion for physicality of the differ- 
entials of the transform. We show in Sect 6 that this physicality close to 
the origin requires a variable c(t) = aE, the normalized cosmic scale factor 
times the normalized Hubble ratio. An infinite number of such coordinate 
systems do exist with these requirements (Sect 7), some of which maintain 
physicality over a large fraction of the universe's space-time. When we use 
a generalized time t we obtain a generalized Equivalency Principle and field 
equation (Sect 8). This allows us to calculate a(t) and thus c(t) (Sect 9). 
The variability of other physical "constants" is discussed in Sect 10. Follow- 
ing the conclusions in Sect 11, we show in the Appendix the derivation of 
some physical coordinates and the results of numerical integration to obtain 
some of them. 

2 The assumptions 

In this paper I believe I can show that a homogeneous and isotropic universe 
requires a physical photon velocity that varies with the cosmic time. Because 
this is such a major departure from the accepted paradigm that the photon 
velocity is constant, I will endeavor to state as clearly as possible what 
assumptions are made to prove this, and then make a stab at additional 
assumptions to handle the major implications resulting from this conclusion. 

Listed below are the four assumptions necessary to derive the variable 
physical photon velocity c(t). The first three are the usual ones for an RW 
universe, but adapted to the possibility of a variable photon velocity; the 
next one determines c(t) by adding the requirement of physicality close to 
the origin of a covariant transform. The assumptions in quotation marks 
are enumerated with italic Roman numerals which will be used throughout 
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the paper when the assumption is being used. The remarks not in quotation 
marks are descriptions, definitions or derivations. 

i. "The universe is isotropic and homogeneous in space." 

The adjective "generalized" will be used to indicate that the cosmic time 
cot is replaced by t, for which dt = c(t)dt. (t has the dimensions of length). 
The generalized RW metric (eq 6) is derivable from i. 

u."We assume the usual interpretation of the RW coordinates applies to 
the generalized RW coordinates." 

The RW metric has time and radial space differentials of dt, a{t)d\ which 
are physical quantities, measurable by clocks and rulers or their technolog- 
ical equivalents, so that the physical photon velocity is c(t) = a(dx/dt) s , 
and the physical velocity of a moving object located at t, x, labeled R, is 
a{dx/dt)R. A galactic point is the terminus of x so that X is a co-moving 
coordinate with which the galactic point stays constant, t is the time on 
a clock co-moving with each galactic point, including the origin at x = 0- 
Each galactic point is in free fall. 

Hi. "The line element ds in the generalized RW metric of ii is an invari- 
ant that can be used to obtain covariant transforms and covariant vectors." 

The alternative line element dr = ds/c(t) could also be used to determine 
the same variable photon velocity, but is less convenient to use in covariant 
relationships. A covariant transform will solve the field equation when RW 
does. 

We will consider a radial world line measured with an isotropic (i) covari- 
ant transform from RW (Hi) whose time and radial coordinates are T(t, %) 
and R(t, x) (eq 7) and whose space origin is attached to the same galactic 
point as x = 0, at which point there is no motion between them, which we 
will therefore call stationary. 

iv. "If the expansion to lowest order in x requires the same variable pho- 
ton velocity for every possible covariantly transformed T and R which have 
effectively 1 a generalized Minkowski metric (eq 21) in this range of Xi an d 
if there exist such transforms, we assume the required photon velocity must 

1 This is given quantitative meaning in Sect 6 
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be the physical c(t) which enters into all of our physical laws throughout 
the universe." 

Essentially, we are assuming that to have physical coordinates the trans- 
formed system must be a stationary inertial system out to some value of x 
with which to measure the expanding galactic points. Expanding to the low- 
est order of \ wm produce the lowest order Lorentz dilation of transformed 
time T, which will exclude both the proper coordinates and a constant c. 
By i, the physical c(t) determined by iv is applicable to the origin placed 
on any galactic point, i.e., everywhere in the universe. 

Additional assumptions are needed to apply this variable physical pho- 
ton velocity to physical laws. I suggest the following: 

v. "The Equivalency Principle be altered to use the generalized Lorentz 
transform instead of the familiar Lorentz transform as the basis for all phys- 
ical laws." 

vi. "Generalized time substituted for normal time makes Einstein's field 
equation applicable to a variable c(t)." 

vii. "The fine structure constant and Rydberg constant (expressed as a 
frequency) remain constant while some other physical "constants" vary with 



3 The generalized RW metric 

Implicit in i is that the concentrated lumps of matter, like stars and galaxies, 
can be averaged to the extent that the universe matter can be considered 
continuous, and that the surroundings of every point in space is isotropic 
and is the same for every point. 

By embedding a maximally symmetric (i.e., isotropic and homogeneous) 
three dimensional sphere, with space dimensions r, 0, and 4>, in a four dimen- 
sion space which includes time t, one can obtain a differential line element 
ds [4, page 403] such that 



ds 2 = g(t)dt 2 - f(t) 



dr 



1 — kr 2 



+ r 2 d6 2 + r 2 sin 6 



(1) 
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a (*) = V /(*) is a cosmic scale factor multiplying the three dimensional 
spatial sphere, so that the differential radial distance is a(t)dx defined by 

The g{t) has normally been taken as g(t) = c 2 = constant, so that c is 
the constant physical velocity of light and t is the physical time on each co- 
moving point of the embedded sphere. In both cases by physical, we mean 
that their value can be measured by physical means like standard clocks 
and rulers, or their technological equivalents. In order to accommodate 
the possibility of c(t) being a function of time, we make g(t) = c(t) 2 . The 
resulting equation for the differential line element becomes a generalized RW 
metric: 

ds 2 = c{t) 2 dt 2 - a{tf[d X 2 + r 2 duj\ (3) 

where 

{sinx, k = 1, 
X, k = 0, (4) 
sinhx, k = -1, 

k is a spatial curvature determinant to indicate a closed, flat, or open uni- 
verse, resp., and dJ 2 = d9 2 + sin 2 9d(p 2 . This is the RW metric generalized 
to a variable photon velocity c(t). 

It will be convenient to introduce the time related quantity t [3], which 
we will call a generalized cosmic time, defined by 

i= f c{t)dt, (5) 
Jo 

so that the line element becomes 

ds 2 = di 2 -a 2 (d X 2 + r 2 duj 2 ). (6) 

4 A covariant transform to measure cosmic dis- 
tances 

The RW coordinates define a differential ruler distance ad\, but the full 
coordinate % ls n °t measured by rulers, but by numbered galaxy positions. 
The proper coordinates T = t, R = ax are an attempt to supply a ruler 
measured distance, but these are not physical in that the resulting metric 
does not become a M metric when expanded from the origin to first order in 
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X {iv). I have tried to find a transformed system which provides coordinates 
that are physical. This is not possible for all values of x (see Sect 7), but is 
possible close to the origin, i.e., out to a small but finite value of %. (There 
are other coordinate systems, such as the full RW system, which describe 
physical phenomena whose coordinates do not use clocks and rulers, but 
my objective was not to find all such systems, but only those that do use 
physical clocks and rulers, which I have thereby called physical). 

The most general line element for a time dependent spherically symmet- 
ric (i.e., isotropic (i)) line element (Weinberg, p335) is 



ds 2 = c 2 A 2 dT 2 - B 2 d/R 2 - IcCdTdR - F 2 (d9 2 + sin 2 9 s d(j) 2 s ) (7) 



where A, B, C, and F are implicit function of T and R, but explicit functions 
of t and x- 

We will look for transformed coordinates which have their origins on the 
same galactic point as x = 0, so R = when x = 0, where there will be no 
motion between them, and where T is t, since the time on clocks attached 
to every galactic point is t (ii), including the origin. Although it doesn't 
directly affect the derivation of the variable photon velocity, in order to have 
a complete transform, we will make 9 S = 9 and 4> s = <fi and make F = ar. 
We will henceforth drop the subscript s. 

We will initially consider only radial world lines so R and T will be 
functions of only t and x : T = T(t, x) and R = R(t, %). So 



where the subscripts indicate partial derivatives with respect to the subscript 
variable. 

We will use as the RW coordinates t, x, 9, 4> and as the transformed 
coordinates T, R, 9, (ft, where dT = cdT, so that the transformed metric 
becomes 



Let us consider a radial point at R on the transformed system. When 
measured from the RW system, it will be moving at a velocity (ii) given by 



dT = T t dt + T x dx, 
dR = R t dt + R x dx, 



(8) 



ds 2 = A 2 df 2 - B 2 dR 2 - 2CdRdf - a 2 r 2 dio 2 , 
and the transform of differentials from (t, x) to (T, R) becomes 

df = T t di + cT x dx, 
dR= \R t di+R x dx, 



(10) 



(9) 




(11) 
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where we will be using V = V/c for convenience. To get the time component 
of the contravariant velocity vector U* = di/ds, we divide eq 6 by dt 2 with 
duo = to obtain 



ds 2 
{ dt' 



«(i) 2 (§) 2 



(12) 



(13) 



To get the radial component, we use the chain rule applied to eqs 11 and 
12: 

dx d\ dt V 
ds di ds 

Since the contravariant velocity vector of a point on the transformed 
system in the transformed coordinates is (1/^4,0,0,0), this will transform 

from the vector = ^ in the RW coordinates (7, ^,0,0) the same as 
df, dR in eq 10: 



1/A = T t7 + ±T xl V, 

= l Rtl + lR xl V. 
Manipulating the second line of eq 14 gives 

aR t 



V 



cR y 



If we invert eq 10, we get 

dt = jj(R x df — cT x dR), 
dx = ±(- l c Rtdf + T t dR), 

where 

D = T t R x - R t T x = T t R x (l + V- 1 x 



aT t 



(14) 



(15) 



(16) 



(17) 



using eq 15. Because of the invariance of ds 2 (Hi), we can enter dt and dx 
of eq 16 into the RW metric (eq 6 ) to find the values of the coefficients of 
eq 9: 

1 - V 2 



A 2 = 



l 

3? 



B 2 = a 



(l + V^) 2 

i-(t) 2 
(i + v c ^f 



and 



C = —. 



T t R 



tn x 



(i + v c ^f 



(18) 
(19) 

(20) 
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5 The physicality criteria. 



In order for T and R to be physical coordinates, we are assuming that the 
general metric of eq 9 near the origin must reduce to the M metric (iv): 

ds 2 = dT 2 -dR 2 -R 2 du 2 , (21) 

I will now try to support this requirement. We have extensive measurements 
on earth that in the absence of gravitation the Lorentz transform between 
inertial systems is accurate and invariant, and keeps the Minkowski metric 
invariant, when c is substantially constant. This suggests that a requirement 
for physicality of an inertial system is that the physical metric to accomo- 
date a variable c(t) must reduce to the Minkowski metric for t practically 
constant. This leads to A = 1, B = 1, C = 0, and R = ar in the transform 
of eq 9. We need the transform to be inertial near the origin because it is 
attached to the inertial galactic point at the origin. We will find in Sect 6 
Eq 33 that the requirements on A, B, and C will automatically yield the 
requirement on R. Thus, for physicality we need the M metric to create a 
stationary inertial system close to the origin (iv). 

In further support of C = 0, note that the radial velocity v p , using ds = 
in the transformed coordinates, is given by 



The two values of v p represented by the ± must be for an incoming and 
an outgoing velocity. Clearly the physical photon velocity cannot be mul- 
tivalued at the same point in space so C must be zero and the metric be 
diagonal, in agreement with eq 21. 
Putting C = in eq 20 we get 

F = -% (23) 

With this, the expressions for A and B in eqs 18 and 19 become: 

A = 1 (24) 
J-t 

B - = (25) 



and v p becomes 



v p = ±A. (26) 
aT t 
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The two different equations for V, eq 15 for distance and eq 23 for time, will 
provide the critical constraint on the expansion of T and R to lowest order 
of x that will determine the variable photon velocity (Sect 6). 

In further support of A = 1 and B = 1, we will use the criteria developed 
by Bernal et al[7] to determine whether two coordinate systems use the same 
clocks and rulers. In their notation for x\ and x 2 translated to the notation 
in this paper, the RW metric has the time measure dt\ and the length 
measure dx\ = a(t)dx and the transformed time measure dt2 = dT and a 
length measure dx2 = dR. Then the criteria for dT physicality become 

dit 2 = T t = d 2 h =t T = j 2 /T t ; (27) 

i.e., A = 1. For dR physicality 

8 lX2 = R x /a(t) = d 2Xl = a{t) XR = J 2 a(t)/R x ; (28) 

i.e., B = 1. Thus, the criteria for physicality of dT and dR lead to the 
requirement 

T t = 1 ,R x = 1 a(t). (29) 

The departure of A and B from unity are measures of how good an approx- 
imation to being physical is a particular diagonal transform (iv). 

The M metric (eq 21) is invariant under an invariant L transformation 
between a frame (dT', dR') moving at a radial velocity iv relative to a second 
frame (dT,dR): 

df = i(t)(df - ^dR), 
dR' = i(t)(-^df + dR), 

where i(t) = 1/^/1 - v 2 /c(t) 2 . Of course, because of spherical symmetry, 
this applies to a frame moving in any direction for differential distances 
measured in the direction of the velocity (i). The M metric and L transform 
become the usual ones for events that occur in times short compared to the 
cosmic time t. Notice that the time dilation (dT'/dT) R = (dT/dT') R > = i 
and distance contraction (dR' / 8R)t = (dRfdR')T' = Y is the same as for 
the usual Lorentz transform. 

6 Transformed physical coordinates require a vari- 
able photon velocity c(t). 

The physicality criteria of iv gives A = 1 = Ti/7 and B = 1 = R x /aj (eq 
29) so that together with the two relations for V eqs 15 and 23 all for small 
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X, we have two relations each for T and R in terms of V. We will proceed 
by using the two relations for R to determine the expansion coefficients of 
R and V, then use the resultant expansion of V in the two relations for T 
to find the requirement for c{t). 

Since V will vanish at the origin (iv), the constant in the power series 
for V is zero: 

-V = Wl (t) X + w 2 (t) X 2 + w 3 (t) X 3 + 0(X 4 )- (31) 
so from R x = 07 (eq 29) we get 

R x = a(l + \v 2 + ^V 4 + ...) = o(l + \w\ X 2 + W2X* + OCX 4 ))- (32) 

If we integrate eq 32 at constant t, noting that R vanishes at x = we 
obtain 

i? = a X + t;Ciw 2 x 3 + -.aw x w 2 X^ + 0(x 5 )- (33) 
b 4 

The first term of eq 33 is the proper distance to which all measurements of 
distance reduce close to the origin [4]. 

Partial differentiation of eq 33 by t at constant x gives 

Rt = ca X + \x*j t {aw 2 ) + -^j t {aw lW2 ) + 0(x 5 )- (34) 

where the dot represents the derivative with respect to t. We can find then 
V from eqs 15, 32, and 34: 

~V = ^ x =«X + f(t)x 3 + 0(X 4 ).-, (35) 

where 

f(t) = -\wla + ^ t {awl) (36) 

By comparison of eq 35 with eq 31, we see that w\ = a, W2 = 0, and 
^3 = f(t). 

We will now use this expression for V to find two relations for Tf. The 
first comes from T t = 7 (eq 29) : 

r t = i + ^ = i + loY + 0(x 4 )-> (37) 

We note that this shows that the galactic clocks t measured at constant 
X run slower than the stationary clocks T as they move away from the 
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origin. When measured at constant R, the stationary clocks run slower, 
tx = 1 + V 2 /2 (eq 16), in accordance with the Lorentz transform. 
We can find an expression for T x , using eqs 23, 37 and 35: 

T x = --T t V = - c [l + id 2 * 2 + 0(X 4 ).-.PX + /(t)x 3 + 0(X 4 )-], (38) 
and multiplying the brackets gives 

T x = - c [*X+ \a 3 X 3 + f(t) X 3 + 0(X 4 )-]- (39) 
By integration with \ a t constant t with T = t at \ = we find 

T = t+ l -^x 2 + 0(x%.. (40) 
If we partially differentiate by t, we get a second expression for T t : 

r < = 1 + ^5<>0<A" (41) 

The equivalence of the two expressions for T t (eqs 37 and 41) out to 
the lowest term in \ {^ v ) requires that there be a variable c(t) given by the 
differential equation: 

° 2 

This differential equation can also be found by using eq 32 with eq 41 to 
calculate v p in eq 26 and setting v p = c(t). 

The reason c(t) appears in this equation arises from T on the stationary 
frame being synchronized by the integration of its differential at constant t. 
Because of a(t), the velocity of the galactic point V varies with time t either 
at constant \ or constant R. The combined constraints on this velocity 
of covariance (—V = cRt/aR x ) and physicality (—V = cT x /aT t , 7 = T t , 
7 = R x /a) requires the photon velocity c{t) also to vary with time t to 
maintain this synchronization. 

To get an explicit expression for c(t), multiply eq 42 by a, change the 
variable dt to da/ca to yield 

^ = ±d{™) (43) 
a ad c 

One can see that c cx d is a solution, so 

c(t) ait) _ . , ,. 

— = = aE, 44 
c ao 
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where a is the normalized scale factor 



a = a/ao, (45) 
and E is the normalized Hubble ratio H(t) 

H 1 a 

E = JT = jr-- (46) 

The subscript denotes the value at the origin at t = to, the present time. 
We will take Co to be unity, so that c{t) will be measured in units of cq. The 
field equation (sect 8) will enable us to evaluate a and E and thus c(t). 

At this point we would like to examine quantitatively how close to the M 
metric our transformed metric needs to be. Let us consider the coefficients 
A, B, and C one at a time departing from their value in the M metric. We 
find for the equation 37 to be valid from the origin out to some small but 
finite x, that \A - 1| << d 2 x 2 / 2 = V 2 /2. B is less critical, \B - 1| << 1, 
since the equality of eqs 37 and 41 depends only on the zero'th order of 
R x in eq 32. For A = 1 and B = 1, eqs 20 and 18 give for small % that 
T t = 1 + \V(V + 2C), requiring \C\ « ax/2 = -V/2. We have also 
assumed that the first term of the expansion is dominant in determining 
c(t), so \c(t) — v p (t)\/c(t) << 1. All of these inequalities should improve 
the closer x g e ts to the origin. However, if we wish to see how far from the 
origin the transform remains physical, we need to examine how large V can 
be and still preserve them. 

The proper coordinates T = t and R = ax move at a velocity V = —ax 
with respect to the galactic points, but without any time dilation. Thus, the 
proper coordinates have A = 1 — (dx) 2 and C = —ax, which do not meet 
the inequality requirements of the previous paragraph. 

Thus, we have shown that if transformed physical coordinates exist with 
a M metric close to their origin, they must have the variable c(t) shown in 
eq 44. We discuss in Sect. 7 and derive in the Appendix an infinite number 
of such transforms. We can therefore conclude by iv that c(t) = aE is the 
physical photon velocity. 

If we had used dr = ds/c(t) as the invariant line element instead of ds 
{Hi) to define the covariant transform close to the origin we would have 
obtained 

dr 2 = dT 2 - -^{dR 2 + R 2 duj 2 ), (47) 
c(ty 

with the same c{t) as above. But the time dependent coefficients of the 
space elements would provide a major complication if we tried to use dr in 
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the field equation, which cannot be simply fixed by incorporating c(t) in dR 
and R. This is a good reason not to use it. 

It should be observed that this variable photon velocity was derived with- 
out the use of the field equation. The existence of the variability is caused 
by the homogeneous and isotropic symmetry assumed for the universe, the 
validity of the covariant transform (including the invariance of the line ele- 
ment ds), and by the concept of physicality extending to the lowest non-zero 
order of \-> but not by the gravitational energy density, whose only effect on 
c(t) is how it affects the rate of change of the cosmic scale factor a. 

7 Transformed physical coordinates exist 

We show in Appendices A, B, and C that there exist an infinite number of 
coordinate systems which satisfy our physicality requirements. Appendix 
A derives the equations for diagonal transforms (C = 0) using physical 
time for all distances (A = 1) and Appendix B using physical distance 
(B = 1). Appendix C describes similarity solutions of these equations for 
a flat universe. Because the equations in A and B are different from each 
other, they show, as we would expect, that it is not possible to have diagonal 
physical R and physical T simultaneously for all values of t, \ (except for 
an empty universe). 

I have been able to find analytic solutions for all the physical time trans- 
forms, but not for all the physical distance transforms. However, because 
there is an interest in having a physical description for distance in the uni- 
verse, we have a special interest in the physical distance transforms. I have 
been able to find a similarity solution for physical distance for the special 
case of flat space with Q = 1. The result is displayed in Figs. 1 and 2. Here 
we have used the field equations with the generalized time (Sect 8) to derive 
the equations for a(t) = (t/to) 1 ^ 2 and c(t) = (toft) 1 / 4 . Fig 1 plots distance 
against the time at the origin (cosmic time t) for galaxies (constant x) an d 
for incoming photons reaching the origin at t/to = 1. The galactic points are 
labeled with their red shift z, determined by the time t of the intersection 
of the photon path with the galactic path z = — 1 + c/a = — 1 + (to / i) 3//4 , 
assuming vii, Notice that light comes monotonically towards the origin from 
all galactic points. This photon path has a slope of cq = 1 close to the origin 
where the distance R and time t are both physical. 

Although the distance uses physical rulers, the coordinate system as a 
whole may not be physical for times shorter than some limit. A reasonable 
limit might be V 2 /2 = 0.25, R = 2.3(t/t ) 3/4 , shown by the dashed line in 
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the figures. The numerical integration gives A = v p /c(t) = 0.953 at this 
physicality limit, so \A — 1| = .047 is much less than V 2 /2. Together with 
B = 1,(7 = for these physical distance coordinates, this shows that the 
assumption that T and R represent a physical inertial coordinate system 
inside this limit is very good. 

Fig 2 plots these distances vs the transform time T. At the emission of 
the photons, T is finite (even for t = 0), presumably the transformed time 
it takes for the galactic point to get out to the point of emission. The slope 
of the light path in Fig. 2 is v p = cq = 1 at T/to = 1, and is only 5% less 
than c(t) = 1.50 at the physicality limit T/to = -40. At the intersection 
of this physicality limit with the photon path that arrives at the origin at 
to, the time t/to = 0.2 and the red shift z = 2.4. Thus, if we have a flat 
universe with £7 = 1, the last 80% of the universe history out to a z of 2.4 
can be treated with physical coordinates T and R. This z is as large as any 
of the Supernovae la whose measurements have suggested an accelerating 
universe. Even out to a z = 7, A - 1 = 0.2 < V 2 /2 = All, B = 1, C = 0, 
Vp/c = A = .8, c(t) = 2, and t/to = -06, so the approximation to M is not 
too bad for a universe 6% of its present age, with a photon velocity at twice 
its present value and its coordinate photon velocity only 20% less than that. 

Actually, the transform is closer to M than keeping only the lowest 
order of the expansion would suggest. Thus, for the plotted physicality 
limit, z = 2.4, the first term of the expansion of 7 — 1 is only 60% of the 
total, while for z = 7, it is only 30%. The locus of points where t = 0, 
labeled horizon, is interesting, but does not represent physical coordinates. 

I also show in Appendix C similarity solutions for physical time for a 
flat universe. These are very similar to the solutions for physical distance, 
but with different numerical limits at t = 0. They are a special case of the 
general solution in Appendix A, but are easier to integrate than the general 
solution. 

8 Revised Equivalency Principle for a variable pho- 
ton velocity. 

To accommodate c(t), it would seem appropriate to revise (v) the Equiva- 
lency Principle as stated by Misner et al[5, page 386] to read "to any and 
every local [generalized] Lorcntz frame, anywhere and any time in the uni- 
verse, all the (nongravitational) laws of physics must take on their familiar 
[generalized] special relativistic forms", viz. L transforms (eq 30) with the 
M metric (eq 21). The use of t or T in the field equation (vi) accomplishes 
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this: 

GV + A<v = ^r^, (48) 

where G^ v is the Einstein contracted curvature tensor determined from the 
coefficients g^ v multiplying the differentials in the metric of eq 6, or of 
eq 9, A is the cosmological "constant" possibly representing some kind of 
vacuum energy density, G is the gravitational "constant" and T^ v is the 
energy-momentum tensor. Locally, this will give the M metric (eq 21) and 
L transform (eq 30) for dT and dR close to the origin, as well as for di and 
adx- 

It will be convenient to use the variables t, r, 9, <f) in the generalized RW 
metric (eq 6): 

ds 2 = di 2 - a 2 {- + r 2 d6 2 + r 2 cos 2 Ode/) 2 ). (49) 

so the g^ v of eq 48 are the coefficients of eq 49, exactly the same as the usual 
RW metric with constant c = 1 in the variables t, r, 9, <fi. We will make the 
usual assumption that the universe is an ideal fluid with an energy density 
of pc 2 and pressure p, so that we can write the two significant field equations 
[5, page 729] for a(t) as 

3d 2 3k 8vrG , . 

— + — - A = — p, 50 

a z a z c z 



and 



da 2 k 8ttG . . 

+2- + — + — - A = —p, 51 

a a z a A c 4 

where the dots represent derivatives with respect to t. 

Now, the first equation can be multiplied by a 3 /3, differentiated, and 
subtracted from da 2 times the second to give 

d Go o. 3G . n x 

Ji^ a) = -^ aap -^ (52) 

where we have allowed for the possibility that G, p, and A may be functions 
of c(t). For small p and A, 

Gpa 3 /c 2 = constant (53) 

If the energy density consists of n particles per unit volume of mass m, so 
p = nm, then the conservation of particles requires na 3 be constant. This 
makes 

Gm 

— 7T- = constant. (54) 
& 
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Even when p and A may not be negligible as in the early universe, it is 
reasonable to assume that for matter density eq 54 will remain valid. 



9 Calculation of a(t) and c(t). 

Following Peebles [6, page 312], we define 

O - fg? <«) 

and 

— k 

= T7272 ( 56 ) 
n a Q 

and 

= (57) 

For very small a there will also be a radiation term whose value will have 
to be reexamined for its dependence on c(t). For this paper it will not be 
included. 

Then the normalized Hubble ratio E in eq 46 is determined by eq 50: 



E = \I^T + ^ + n A . (58) 
V or cr 

which allows us to evaluate c(t) = aE. The S7s are defined so that 

n + n r + n A = 1. (59) 

At t = to: a = 1, E = 1, and co = 1. 
The cosmic time t becomes 

f a da l' a da 

c ° H ° t= l 5s=y„ sf- (6o » 

where Co = 1 is written explicitly for clarity. For £7 = 1: collet = a 2 /2, 
whence c H t = 1/2, a = (t/t ) 1/2 , and c = a' 1 / 2 = (to/t) 1 / 4 . 

For experiments attempting to measure the variation of the photon ve- 
locity at the origin at the present time, the derivative of c(t) (eq 44) will be 
more useful: 



1 



cqH 



i del 3 n 
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Notice that this fraction is negative when matter dominates, and goes from 
zero at zero density to —1/2 at the critical universe density. A vacuum 
energy density opposes the gravitational effect of matter, and as the vacuum 
energy density increases, the slope goes to zero before reversing from a 
decreasing function of time to an increasing function. 



10 Other physical constants. 

Many physical quantities depend on c. The fine structure constant aj in SI 
units [8] is 

a > = 4^ft? < 62 > 

and the Rydberg constant expressed as a frequency is 

2 m e c 2 e 4 m e 

R - c = a f^n=4^W (63) 

Note that the 4-ireo is often omitted in the fine structure constant since it is 
unity in Gaussian coordinates, but it is essential here if we are to consider 
a variable c{t). 

Atomic clocks have reached a stability and precision to demonstrate that 
a/ and RooC remain essentially constant with time [9] [10] (vii). From this 
we can conclude that 

^ oc c(t), (64) 

and 

me 1 fan 

T a # (65) 

Another relation we can count on comes from electromagnetic theory: 

eo^o = -77T2, (66) 

c{ty 

(with c{t) not normalized here). 

Lastly, if the use of generalized cosmic time in the field theory is correct, 
we obtain (eq54): 

Gm oc c{tf. (67) 

If e and H are constant, we could conclude that eo = l/c(t), fiQ = l/c(t), 
m oc l/c(t) 2 , and G oc c(t) , which would keep constant the mass energy 
mc 2 and gravitational energy Gmim2/2r 12 , but this is only one possibility. 
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The actual variation will depend on comparison with experimental data, 
such as supernovae la, cosmic background radiation, gravitational lensing, 
and dynamical estimates of galactic cluster masses. For this, of course, we 
need a theory of the c(t) dependence of each of these. 

11 Conclusions 

I have investigated what coordinate systems, other than the Robertson- 
Walker (RW) system, there might be to describe a homogeneous and isotropic 
universe that are physical close to the origin, i.e., are derived by a covari- 
ant transform from the RW coordinates and use physical clocks and rulers 
with a Minkowski metric. I have shown that all such systems must have the 
same variable photon velocity close to their origin, which is consistent with 
the RW system also having the same variable photon velocity. Since this 
is determined by physical coordinates, it is reasonable to suppose that this 
represents the physical photon velocity. I find that the analytic expression 
for it is c(t) = aH, the cosmic scale factor times the generalized Hubble 
ratio in appropriate units. For a homogeneous universe, the origin can be 
placed at any point in the universe, so this variable photon velocity should 
be applicable to all physical processes everywhere. 

The existence of the variability is caused by the homogeneous and isotropic 
symmetry assumed for the universe, the validity of the covariant transform, 
and by the concept of physicality of the transformed coordinates close to 
the origin. If we deny the variability of the photon velocity, we must find a 
flaw in one of these concepts. It is not caused by the gravitational energy, 
whose only effect on c(t) is how it affects the rate of change of the expansion 
coefficient a. 

A variable c(t) no longer keeps the local Lorentz transform invariant. 
However, we can make it invariant by using a generalized time for which the 
differential is c(t)dt. We can then create a covariant derivative and a local 
generalized Minkowski metric in order to use Riemann algebra in Einstein's 
field equation. With this we find that the Hubble ratio defined with this 
generalized time has the same dependence on the cosmic scale factor a as it 
does for a constant photon velocity. With the use of this generalized time 
in the field equation, I have found transforms that remain physical over a 
large fraction of the universe space-time 

Some other physical "constants" must also be changing with cosmic 
time, but its hard to say which ones. I have suggested some constraints on 
this variability, including keeping constant the fine structure and Rydberg 
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constants. One possibility might be to keep constant the electron charge, 
Planck's constant, mass energy, and gravitation energy, with vacuum per- 
mittivity and permeability each varying as l/c(t), but other variations are 
also possible. 



Appendix A: Transformed diagonal coordinates with 
physical time 

For diagonal coordinates with physical time at all t and %, A = 1. Thus, eq 
24 becomes 

T t = 7- (68) 

Combined with C = 0, this automatically guarantees the Lorentz time dila- 
tion (eq 27). We need only find a transform using it for which B — > 1 close 
to the origin. 

We proceed by finding a differential equation with only V as the de- 
pendent variable. Thus, we can write a formula for T, using this and eq 
23: 

T = t+ f X T x d X t = t + [ X (-- 7 V)dxt, (69) 
Jo Jo c 

where we have used the boundary condition that at x = 0, T = t, and the 

symbol dxt signifies integration with x at constant t. It can be partially 

differentiated with respect to t (giving 7) and then with respect to x an d 

with the use of eq 11, noting that = ^VdV and 1 + V 2 ^ 2 = 7 2 , we 

obtain a PDE for V: 

This can be rewritten as 

dV R d^) R 



V{l-V 2 ) 



(71) 



where the subscript on the partial differential indicates the variable to be 
held constant. This can be integrated with an integration constant hiK. 
Since the integration is done at constant R, then k = n(R), and inversely, 
R = R(k). Integrating eq 71, we get 

V = - " ■ (72) 



4 + K 2 
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where the sign of k will be positive for an expanding universe, where the x 
points will stream out radially past a point at R. 

At this point, R is an unknown function of k. The various possible 
coordinate systems which solve our PDEs are characterized, in large part, 
by the function R(k). But for all, in order for V to vanish when R = 0, k 
must also; so always 

k(0) = 0. (73) 

We note that as long as k(R) remains finite, V goes to —1, and V goes to 
—c(t), for a(t) = 0, i.e. for t = 0, the horizon. 

Let us now look at lines of constant k(R), i.e. constant R, in t, x space. 
Eq 11 can be integrated for x with use of eq 72 at constant k to give the 
following: 

X(t,n)=r CK fl K ■ (74) 

For an open universe, we have set the upper limit at oo, because we expect 
that if R is kept constant the galactic point x that will be passing any given 
R will eventually approach zero as RW time t approaches infinity. This 
limit may be valid even for a closed universe, since c(t) goes to zero for long 
times, making the time go to infinity perhaps before the first maximum of 
a periodic R. 

At this point, we have obtained V = V*(t,n) from eq 72 and have also 
obtained the function x(t> K )- We can in principle invert eq 74 to obtain 
k in terms of t and x : K = K(t,x)- This gives us the velocity function 
V(t,x) = V*(t,K(t,x)- If the function R(k) were known, we would then 
also have R(t,x) = R{K(t, X ))- 

T(t, x) can be found by noting from eqs 23 that 

aV m aV 
T x = -—T t = -— 7 = «. (75) 

By substituting eq 75 into eq 69, and integrating over k instead of x by 
dividing the integrand by the partial of eq 74 with respect to k, we find an 
expression for T(t, x): 



T(t ' X) =t + It 



l + ^# 



ds K , (76) 



where k is put equal to K(t, x) after integration at constant k in order to 
get T(t, X ). 
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This completes the solution. Since k(R) can be any function that van- 
ishes at the origin, there thus exist an infinite number of solutions for our 
transformed coordinates with A = 1, C = 0. 

To determine physicality, we will next find close to the origin. 

R x = k x /k'(R) can be written in an inverted form by taking the derivative 
of eq 74 with respect to k at constant t: 



Ry 



= K 



dK 



tJ 



K 



dn 



AR) J 



c 2 ds K 



t a 2 (l + ^1)3/2 



(77) 



To be physical B = ^a/R x — > 1 as R approaches 0. Putting 7=1 and 
k(0) = in eq 77, multiplying by a, setting the resultant equal to one, and 
changing the integration variable from t to a(t) gives 



1 

a 



K '(0) 



cda 



(78) 



remembering that the dot indicates differentiation by t. k'(O) is a constant 
to be determined by c(to) = 1. Note that the integral of eq 78 is independent 
of the functional form of n(R), and is therefor the same for all k(R). In fact 
it is an integral equation for c(t). By differentiation of both sides of eq 78 
by a, we can obtain 

cit) = ^a, (79) 

which, as we should expect, is the same c(t) (see eq 44) we showed for all 
physical coordinate systems for k'(0) = do- 

Although there is a physical transform for every k(R), we would like to 
find the one that extends the physicality limit out from the origin as far as 
possible. To do this we will try to find a n(R) which eliminates the first 
derivative of B close to the origin. From eqs 25 and 77 we see that 



B(t,R) = ^-=-ya K '(R) f 
K x Jt 



c 2 ds K 



t a 2 (l + ^1)3/2 



= W[R)J t - 



c 2 ds K 
a 2 7 3 



(80) 



where we have used the relation 7 = (1 + ^§^) 1//2 . Since 
we can differentiate 

dB 



87 
dR 



J t 



(t,0) = 0, 



dR 



(t 0) = ^ 



(81) 



By setting k"(0) = 0, this eliminates the first derivative of B(t, 0) to give 



k = k'(0)R 



(82) 
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close to the origin. Any function k(R) that is linear with R close to the 
origin will extend the physicality limit the farthest away from the origin. 
The simplest, of course, is to keep k linear for all R. 

It is interesting to consider the limiting case of a zero density universe: 
Q = 0,Q r = 1, (IqHo = 1 (eq 56). Eq 58 makes a = 22 for all t, %. 
Integrating gives a = t. Eq 44 makes c = 1. Eq 74 gives % = csch _1 t/K, or 
k = K(t, x) = t sinhx- We can then find from eq 72 that V(t, x) = ~ tanhx 
and from eq 12 that 7 = cosh \ so that 

R x dRK x dR , . 

c= =1 f = l- = 1 - ^ 

7a an 71 an 

Thus the physicality condition is met for all R with R = K and A = 1,2?= 1, 
so that the complete transform with eq 76 becomes 

R = tsinhx, ^ 
T = t cosh x- 

These coordinates have been known ever since Robertson [11] showed that 
this transformation from the RW co-moving coordinates at zero density 
obeyed the Minkowski metric. What is new is that this solution was derived 
from the equations we obtained for our physical time transforms with A = 1. 
It can also be obtained from the physical distance transforms (B = 1) since 
eqs 70 and 89 for V become identical with V t = and a/c = a = t. It is the 
only known rigid physical coordinate system for all times and distances in 
a homogeneous and isotropic universe. 

Appendix B: Transformed diagonal coordinates with 
physical distance 

For diagonal coordinates with physical distance for all t and x> B = I, that 
is 

R x = 07 (85) 

By integration we find 

R = af X 1 d X u (86) 
Jo 

and partial differentiation with respect to t gives 



rx rx 
R t = ca / -fdxt + a ifixt- (87) 
Jo Jo 
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We can then find V from eq 15 as 



07 



1 

07 



ca 

Jo 



x rx 

idxt + a j t dxt 
Jo 



(88) 



This is an integral equation for V . It can be converted into a partial differ- 
ential equation by multiplying both sides by 7 and partial differentiating by 



7 



v x + -vv t 

c 



= —a = 



1 da 
c dt 



(89) 



Note that this is substantially different from the eq 70 for V that we ob- 
tained for physical time This means that it is not possible to find diagonal 
transforms with both physical time and physical distance for all values of 
t and x (except for $7 = 0). It is possible to have either one or other be 
physical at all t and x with the other being physical only close to the origin. 



Appendix C: Similarity solutions for flat universe, 

tt = 1 

I have not found it possible to obtain explicit general integrations of eq 89 as 
we did for physical time. However, I have found similarity integrations for 
the special case of = 1 where a = a (tAo) 1/2 and c = (t A) 1/4 - To simplify 
notation let us normalize time to t/to — ► t, a/a^ — > a, and x°o/co^o — * mx, 
T/to -» T, R/c t -» R, and let W = —V. Eq 89 then becomes 

Lw x - t^WWt = V 1/4 (l - W 2 ) (90) 

This can be converted into an ordinary differential equation (ODE) by let- 
ting 

u = x/t 1/4 (91) 

so that eq 90 becomes 

W'(± + ^f) = \(l-W% (92) 

where the prime denotes differentiation by u. 

Similarly we can find ODE's for T and R by defining: 

T = tq(u), (93) 
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and 

R = mt 3/4 s(M), (94) 
where q(u) and s(u), from eqs 23 and 15, are given by the coupled ODE's: 

and 

m 4 4 

It is useful to find that s = j(u+4W/m)/3 and ^ = j/T t = j(l+muW/4)/q. 

For small values of u, W = mu/2, q = 1 + m 2 u 2 /4, s = u, and i? = 
^i/2 x _ ax _ -^y e a j go no t e that T t = 1 + W 2 /2, confirming that these 
coordinates have physical time close to the origin, justifying c(t) = i -1 / 4 . 
An alternate approach would be to start with c(t) unknown, but of the 
form c = t~ b . This will also yield a similarity integration of eq 89 with 
T t = 1 + 2bW 2 for small u. In order to make the time physical close to the 
origin, Tt has to equal 1 + W 2 /2 so that b = 1/4, confirming that for these 
diagonal coordinates to be physical at the origin c 2 oc da/dt, or c oc da/dt. 

For m = 1 as t -> 0, u -> oo, and W 2 -> 1 - 760/(n + 4) 4 , g -» 
.00132(n + 4) 4 , and s -> .0121 (u + 4) 3 . T and i? both remain finite at this 
limit with T -> .00132x 4 , and ,R -> .0121x 3 . The fact that T does not 
go to zero when t goes to zero results from synchronizing T with t at the 
origin, so that at finite R, where clocks run slower than t, they have different 
starting points. The distance R and time T can be found from the numerical 
integration of the coupled ODE's. The paths of galactic points are those for 
constant x. The path photons have taken reaching the origin at t\ is found 
by calculating x p vs t and using the transform to R. Thus, for Q = 1 

x - [ t °dt = ±(t\ /A -t 1 '*) (97) 

Jt! a 171 



Galactic and photon paths for m = 1 are shown in Fig. 1 and 2. An 
approximate upper limit of physicality is shown by the heavy dotted line: 
W 2 /2 = 0.253, u = 2.0, R = 2.30t 3 / 4 = 1.35T 3 / 4 . R vs T at t = provides 

3 /4 

a non-physical horizon: R^ = 1.747T, ' . 

There are also similarity solutions for physical time. With the same 
simplifications as above, the ODE for W is 

W'{- + j) = - A W{l-W 2 ), (98) 
m 4 4 
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with the ODE's for T.R.q, s, x p ) the same. These all are the same as physical 
distance for small u, but differ numerically at large u. Useful relations 
for physical time are obtained from the general solution in Appendix A: 
R = 2k/ 'm, q = 7(1 + muW/4), s = 2jW/m, and B = j(u + 4W/m)/3s = 
2(l + um/4W)/3. 

There are an infinite number of similarity solutions obtained by different 
ratios of x to x, which are still physical close to the origin. They all give 
finite limits of T and R at t = 0, but with different numerical values. I have 
numerically integrated several of both kinds of these without significant new 
insights since they all behave the same at the same rau for small u. 
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Figure 1: Physical distances for $7 = 1 
Distance R/c^tQ plotted against the normalized time on clocks at the origin 
(t/to) for various galaxy paths (labeled by their red shift z) and for the light 
path which the photons take after emission by any galaxy that arrives at the 
origin at t/to = 1. Notice that the slope of this light path close to the origin 
is Co = 1, where the coordinates R and t are both physical. The dashed line 
in Fig. 1 shows the approximate upper limit of physicality, where both R 
and and the transformed time T are physical. 
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Figure 2: Physical distance vs transformed time for = 1 
Distance R/coto plotted against the transformed time (T/to) on clocks at- 
tached at R for various galaxy paths (labled by their red shift z) and for the 
light path that arrives at the origin at T = to. The horizon is the locus of 
points where t = 0. The dashed line in Fig. 1 shows the approximate upper 
limit of physicality for the transformed coordinates. The slope of the light 
path is very close to c(t) out to the limit of physicality. 
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